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1.  Introduction 


This  is  a  first  report  on  diffusion  experiments  done  in  summer,  1993,  at  the  EPA 
Fluid  Modeling  Facility  (FMF)  by  Mark  Piper  with  the  generous  guidance  and  assis¬ 
tance  of  their  staff,  including  FMF  director  Bill  Snyder  and  FMF  scientist  Bob  Lawson. 

These  experiments  are  part  of  a  project  whose  objectives  are  to  improve  our  un¬ 
derstanding  of  diffusion  in  convective  turbulence  and  to  develop  better  models  of  this 
diffusion  for  a  wide  range  of  applications.  Toward  those  ends,  two  diffusion  problems  axe 
being  studied  in  the  FMF  convection  tank.  The  first  concerns  “transport  asymmetry”, 
the  difference  in  the  diffusion  properties  of  scalars  introduced  at  the  bottom  and  top  of 
convective  boundary  layers.  The  second  area  of  study,  the  diffusion  of  highly  buoyant 
plumes  in  the  convective  boundary  layer,  will  be  carried  out  under  subcontract  by  Dr.  J. 
C.  Weil  of  CIRES  at  the  University  of  Colorado. 


2.  Theory  of  Top-Down,  Bottom-Up  Diffusion 

During  the  1980s  the  numerical  technique  called  large-eddy  simulation  (LES)  pre¬ 
dicted  striking  differences  between  the  “top-down”  and  “bottom-up”  vertical  diffusion 
from  very  large  area  sources  within  the  convective  boundary  layer  (CBL).  Top-down  dif¬ 
fusion  is  defined  as  that  process  in  which  the  vertical  turbulent  flux  of  the  scalar  is  zero 
at  the  surface  and  maximum  at  the  mixed-layer  top;  in  bottom-up  diffusion  the  flux  is 
zero  at  the  top  and  maximum  at  the  surface.  For  example,  ozone  introduced  into  the 
boundary  layer  by  turbulent  entrainment  of  the  overlying  atmosphere  undergoes  top- 
down  diffusion.  The  bottom-up  process — a  scalar  flux  at  the  surface  but  not  at  the 
mixed-layer  top — probably  does  not  exist  in  isolation  from  the  top-down  one.  This  is  be¬ 
cause  the  entrainment  process  in  conjunction  with  the  change  in  the  scalar  concentration 
across  the  inversion  “reflects”  the  bottom-up  diffusion,  generating  a  flux  at  the  top  as 
well.  For  example,  solar  heating  generates  a  positive  temperature  flux  at  the  surface,  but 
the  turbulent  entrainment  of  the  capping  inversion  brings  down  warmer  air  and,  hence, 
generates  a  negative  temperature  flux  at  the  mixed-layer  top.  Thus,  a  heated  surface 
drives  a  combination  of  bottom-up  and  top-down  diffusion  of  temperature,  and  an  evap¬ 
orating  surface  drives  a  combination  of  bottom-up  and  top-down  diffusion  of  water  va¬ 
por. 

It  is  the  linearity  of  the  scalar  diffusion  process  that  makes  the  top-down,  bottom- 
up  concept  a  useful  one.  To  illustrate,  let  us  represent  a  scalar  field  c  in  the  CBL  as  c.  — 
C  +  c,  the  sum  of  ensemble-mean  and  fluctuating  parts.  One  can  decompose  each  of 
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C  -  C  \  0 


these,  in  turn,  into  top-down  and  bottom-up  components: 

C  =  C*  +  Cb',  c  =  C(  ■+•  Cj,. 

The  mean  gradient  and  turbulent  flux  are  thus  the  sum  of  top-down  and  bottom-up 
components: 


dCt  dCb  _ 

“3—  +  -X—;  CIC  : 

dz  dz 


CtW  +  CbW. 


Wyngaard  and  Brost  (1984)  and  Moeng  and  Wyngaard  (1984)  isolated  the  statistical 
properties  of  top-down  and  bottom-up  diffusion  from  LES  data.  They  found  that  the 
top-down  eddy  diffusivity 

*  -  aS  (3) 

is  well-behaved,  having  a  maximum  magnitude  of  about  0.1iu*z,-  in  midlayer.  Here  w»  is 
the  convective  velocity  scale  (fiQoZi)1/3,  where  /?  is  the  Boussinesq  buoyancy  parameter, 
Qo  is  the  surface  temperature  flux,  and  z*  is  the  boundary-layer  depth.  The  bottom-up 
diffusivity  is  everywhere  larger  in  magnitude  and  has  a  midlayer  singularity.  This  means 
that  the  well-documented  failure  of  eddy  diffusivity  closure  for  scalars  in  the  convective 
boundary  layer  (e.g.,  Deardorff,  1966)  is  confined  to  the  bottom-up  process. 

« 

Transport  asymmetry  has  yet  to  be  confirmed  through  direct  observations  in  the 
atmosphere.  It  appears  that  a  specially  designed  field  program,  with  scatter  minimized 
through  optimally  chosen  scalars  and  long  flight  legs,  is  required  to  test  the  LES  results 
in  the  atmosphere.  One  of  our  objectives  is  to  test  the  hypothesis  first  in  the  laboratory. 


3.  Operation  of  the  FMF  Convection  Tank 

Fluid  modeling  allows  one  to  simulate  certain  geophysical  flows  in  the  laboratory 
without  precisely  matching  the  values  of  all  dimensionless  parameters  (Snyder,  1981). 
Deardorff  and  Willis  began  a  long  series  of  studies  in  a  laboratory  convection  tank  1  m 
on  a  side,  using  water  as  the  fluid,  in  the  late  1960s;  they  soon  perceived  it  as  a  model 
of  the  unstable  planetary  boundary  layer  (Willis  and  Deardorff,  1974).  Research  find¬ 
ings  from  this  tank- continued  for  nearly  two  decades  (Deardorff  and  Willis,  1985).  Their 
early  diffusion  studies,  particularly  after  they  were  corroborated  through  LES,  revolu¬ 
tionized  the  way  we  perceive  convective  plume  dispersion  (Lamb,  1982).  These  labora¬ 
tory  experiments  have  provided  much  of  the  present  basis  for  short  range  atmospheric 
dispersion  modeling  under  convective  conditions  (Weil,  1988). 

The  original  Deardorff- Willis  convection  tank  is  now  located  in  the  Fluid  Modeling 
Facility  of  the  Atmospheric  Research  and  Exposure  Assessment  Laboratory.  The  tank 
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sidewalls  are  constructed  of  1.9  cm  thick  plexiglass,  and  the  bottom  is  a  1.27  cm  thick 
anodized  aluminum  plate.  The  horizontal  inside  dimensions  of  the  tank  are  122  cm  x 
124  cm;  the  depth  is  45  cm.  In  the  experiments  done  in  the  summer  of  1993  we  filled  the 
tank  to  35  cm  depth  with  filtered,  deaerated  water.  To  minimize  heat  loss,  5  cm  slabs  of 
styrofoam  are  attached  to  the  sidewalls  of  the  tank  and  smaller  pieces  of  2  cm  styrofoam 
are  floated  on  the  water  free  surface.  To  heat  the  lower  surface  of  the  tank,  144  electric 
heaters  are  affixed  to  the  underside  of  the  aluminum  plate  at  the  tank  bottom.  They  are 
constructed  of  a  metal  foil  insulated  with  rubber  and  coated  with  a  pressure-sensitive 
adhesive  on  one  side.  Up  to  6  kW  total  power  can  be  dissipated  through  these  heaters. 

To  simulate  the  inversion  layer  that  lies  above  the  daytime  CBL,  a  stably  stratified 
temperature  profile  was  formed  by  passing  a  grid  of  heating  wires  through  the  upper  20 
cm  of  water  in  the  tank.  The  heating  grid  consists  of  four  102  m-long  circuits  of  2  mm 
teflon-coated  wire,  each  circuit  having  a  computed  resistance  of  19  ohms.  A  current  of 
2.9  amps  was  applied  to  the  wires,  giving  a  total  power  output  of  over  600  W. 

In  our  experiments  we  heated  the  water  between  15  and  35  cm  to  create  an  inver¬ 
sion  of  strength  0.5  C  per  cm.  The  depth  of  this  stratified  layer  satisfies  the  criterion  set 
by  Deardorff  and  Willis  (1985)  that  the  aspect  ratio  of  tank  sidewall  separation  to  layer 
depth  be  at  least  5:1.  We  used  two  scalars,  temperature  and  food  dye,  in  our  diffusion 
experiments. 

Measurements  of  temperature  diffusion 

We  measured  horizontally  averaged  mean  temperature  with  a  resistance-wire  grid 
formed  by  stringing  33  meters  of  0.762  mm  teflon-coated  alumel  wire  across  an  alumi¬ 
nium  frame.  The  grid  frame  was  constructed  from  1”  x  1/8”  aluminum  strips  and  mea¬ 
sures  120  cm  on  a  side.  We  obtained  vertical  profiles  of  mean  temperature  by  traversing 
the  grid  vertically  through  the  tank  at  an  average  speed  of  0.64  cm  s-1.  A  potentiome¬ 
ter  marked  the  vertical  position  of  the  grid  in  the  tank.  The  grid  position  in  centimeters 
was  interpolated  from  a  two-point  calibration  table. 

We  used  the  voltage  deflection  across  a  Wheatstone  bridge  to  measure  the  RTD  grid 
wire  resistance,  which  varies  approximately  linearly  with  temperature.  A  calibration  ta¬ 
ble  was  compiled  by  measuring  the  bridge  output  voltage  against  a  thermistor  reading; 
temperature  was  then  interpolated  from  this  calibration  table.  A  10  Hz  lowpass  filter 
was  used  to  remove  noise  from  the  output  voltage  signal. 

A  nominal  value  of  3  kW  of  power  was  dissipated  through  the  heaters  at  the  lower 
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surface  of  the  tank.  This  rate  of  surface  heating  provides  a  turbulent  Reynolds  number 
based  on  layer  depth  2,  and  velocity  scale  w,  of  at  least  1000. 

During  a  data  retrieval  run,  output  voltage  was  sampled  continuously  at  56  Hz.  A 
14-point  block  average  was  immediately  applied  to  the  raw  data,  giving  an  effective  sam¬ 
pling  frequency  of  4  Hz.  Data  retrieval  runs  varied  in  length  from  10  to  50  min.  Elapsed 
time,  grid  position,  bridge  voltage  and  temperature  data  were  written  to  a  ASCII  file  at 
a  rate  of  4  records  per  second. 

Measurements  of  dye  diffusion 

We  introduced  the  top-down  scalar,  food  dye,  as  follows.  The  tank  was  initially 
filled  to  a  depth  of  32  cm  with  filtered,  deaerated  water  and  allowed  to  cool  to  nearly 
21  C.  We  mixed  400  ml  of  dye  with  water  for  a  nominal  45  liters  of  solution.  This  solu¬ 
tion  was  kept  at  40  C.  To  induce  thorough  mixing,  the  dye  solution  was  introduced  into 
the  tank  at  a  point  slightly  below  the  water  free  surface  at  an  average  flow  rate  of  25  ml 
s_1.  While  the  dye  was  being  injected  into  the  tank,  the  stratification  system  was  also  in 
operation,  aiding  in  the  mixing  of  the  dye  solution. 

After  the  solution  was  completely  siphoned  into  the  tank,  the  dye  layer  extended 
from  15  cm  to  35  cm,  coinciding  with  the  layer  of  stable  temperature  stratification  in 
the  tank.  A  small  amount  of  dye  tended  to  leak  through  the  stable  inversion  layer;  be¬ 
ing  negatively  buoyant,  it  dropped  to  the  surface  of  the  tank.  These  dye  streamers  were 
quickly  mixed  and  diluted  shortly  after  the  surface  heating  of  the  tank  was  initiated. 

We  sampled  dye  concentration  by  siphoning  fluid  through  plastic  tubing  at  16  dis¬ 
crete  levels  at  five  locations  in  the  convection  tank.  The  fluid  was  collected  in  test  tubes 
and  analyzed  with  a  photocolorimeter.  The  measured  transmissivity  of  a  light  beam 
through  the  sampled  fluid  was  compared  against  a  calibration  table  of  measured  stan¬ 
dards.  This  system  is  appealing  because  of  its  simplicity  and  because  the  sampling  rakes 
fit  between  the  RTD  grid  wires,  allowing  temperature  and  dye  concentration  to  be  sam¬ 
pled  simultaneously. 


4.  Top-down,  Bottom-up  Diffusion — Preliminary  Results 

As  we  discussed,  a  scalar  introduced  at  the  surface  undergoes  both  bottom-up  and 
top-down  diffusion,  the  latter  being  induced  by  the  entrainment  process  in  the  capping 
inversion  that  marks  the  top  of  the  convective  layer.  For  that  reason  there  seems  to  be 
no  way  to  produce  “pure”  bottom  up  diffusion  experimentally.  Because  the  bottom  sur- 
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face  can  be  inpenetrable  to  some  scalars  (e.g.,  density  of  a  diffusing  constituent)  “pure” 
top-down  diffusion  experiments  can  be  done,  however.  In  principle,  using  two  different 
scalars,  either  simultaneously  in  one  flow  or  sequentially  in  two  statistically  identical 
flows,  allows  one  to  deduce  the  mean  transport  properties  of  the  top-down  and  bottom- 
up  processes. 

The  mean  temperature  T(z,t)  in  the  convecting  fluid  layer  evolves  according  to 


dT  ddw 
dt  dz 


=  0. 


(4) 


Figure  1  shows  the  measured  mean  temperature  versus  height  and  time.  From  these 
data  we  calculated  the  vertical  profiles  of  dT/dt  shown  in  Figure  2.  From  (4)  this  gives 
the  divergence  of  the  turbulent  temperature  flux.  Integrating  this  flux  divergence  over  z 
yields,  from  (4)  and  the  lower  boundary  condition, 


_  _  fz  dT  _  _ 

8w(z)  =  8w0-J  —dz'~(l-z/zi)8w0  +  (z/zi)8wi,  (5) 

since  the  profile  is  known  to  be  essentially  linear  over  the  mixed  layer.  Since  tempera¬ 
ture  undergoes  both  top-down  and  bottom-up  diffusion,  both  the  surface  flux  w80  and 
the  top  flux  w8i  are  nonzero. 

In  principle  the  surface  flux  8w0  =  Q0  is  known  from  the  thermal  energy  balance 
for  the  lower  plate  of  the  convection  tank,  so  the  top  flux  8wx  can  be  found  from  (5). 

For  reasons  that  axe  not  yet  clear,  we  found  that  this  gave  unrealistic  flux  profiles  near 
the  fluid  layer  top.  Instead,  we  allowed  the  surface  flux  to  be  determined  by  (5)  under 
the  condition  that  the  flux  in  the  undisturbed  flow  deep  in  the  inversion  layer  was  zero. 
This  gave  the  flux  profiles  in  Figure  3.  The  height  where  the  flux  reaches  its  negative 
naximum  is  taken  as  z,-;  it  is  indicated  on  each  of  the  profiles  in  Figure  3. 

Figure  4  shows  on  an  expanded  temperature  scale  a  typical  temperature  record  from 
the  traversing  probe.  Since  this  record  is  from  a  moving  probe,  time  varies  with  position 
along  the  profile.  Converting  this  record  to  a  single  time  yields  the  temperature  profile 
in  Figure  5.  The  noise  level  in  this  profile  is  encouragingly  low.  The  resulting  vertical 
gradient  of  temperature,  scaled  with  z,-  and  the  surface  temperature  scale  T»  =  Q$/wm,  is 
shown  in  Figure  6.  Also  plotted  there  is  a  profile  from  large-eddy  simulation  of  a  convec¬ 
tive  boundary  layer.  The  agreement  is  good  except  near  the  surface.  We  are  currently 
investigating  the  reasons  for  the  discrepancy  there. 

The  mean  concentration  Ct  of  the  food  dye,  entrained  from  the  capping  inversion 
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layer  by  convective  eddies  in  the  growing  convective  layer  below,  evolves  as 

dct  dm  _ 

dt  dz 

LES  experiments  (Moeng  and  Wyngaard,  1984)  indicate  that  during  this  process  the 
scalar  flux  profile  is  nearly  linear,  so  that 


ctw(z)  =  Ctw(Zi )  =  ctwu 


and  the  top  flux  ctWi  does  not  depend  significantly  on  time.  Thus,  the  mean  gradient  is 
quasi-steady, 

Figure  7  shows  the  measured  mean  concentration  Ct(z,  t)  at  16  heights  z  and  five 
times  t.  The  calculated  time  derivative  dCt/dt  is  shown  in  Figure  8.  The  deduced  top- 
down  flux  divergence  is 

dm  dCt  ,Q, 

dz  dt  ’  {  ' 

and  since  by  definition  ctw(0)  =  0,  we  can  integrate  this  to  get 


ctw(z)  =  {z/zi)  ctw1  =  -  J  -gf-dz'  =  -t 


The  resulting  flux  profiles  at  the  five  times  are  shown  in  Figure  9. 

Figure  10  shows  the  measured  vertical  gradient  of  mean  dye  concentration  nondi- 
mensionalized  with  z,-  and  the  concentration  scale  c*  =  cTtZ^/ia*.  Since  the  dye  diffusion 
is  a  top-down  process,  this  nondimensionalized  dye  gradient  is  precisely  the  top-down 
gradient  function  gt: 

/  /  \  dCtZiW *  ,,,N 

9t(z/zi)  =  — ■£-=-  ■  (11) 

uZ  CfW i 

The  LES  result  for  gt  (Moeng  and  Wyngaard,  1984)  is  shown  in  Figure  10  along  with 
the  dye  data.  The  agreement  is  fairly  good. 

We  estimated  the  bottom-up  gradient  function  as  follows.  By  the  linearity  of  diffu¬ 
sion  and  the  similarity  hypothesis  we  have 

dT  dTb  ,  dTt  0wo  0wi 

TP  —  - 1-  ~x—  — - 9b - gt-  (12) 

UZ  OZ  OZ  WmZi  W+Zi 
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Having  the  top-down  gradient  function  gt  from  the  food-dye  results,  we  used  this  equa¬ 
tion  to  infer  the  bottom-up  function  gt  from  the  temperature  data: 

w*Zi  ( dT  9w\  \ 

9b  —  — ■==-  (  H - 9t  )  ■ 

6w0  \oz  w*Zi  ) 

Figure  11  shows  the  resulting  gt  estimates  and  the  corresponding  LES  results.  The 
agreement  is  somewhat  poorer  that  for  gt,  Figure  10,  principally  because  of  the  system¬ 
atically  low  measured  gradient  values  near  the  surface.  Nonetheless,  the  LES  results  of 
transport  asymmetry — that  gt  at  z/z(  is  systematically  larger  than  gt  at  (1  —  z/zf),  so 
that  the  bottom  up  diffusivity  Kt  exceeds  the  top-down  one  Kt — is  reproduced  by  the 
tank  experiments. 
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Figure  3.  Deduced  vertical  temperature  flux  profiles  for  three  traverses. 
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Figure  4.  An  expanded  record  of  mean  temperature  versus  height  and  time  for  traverse 
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Figure  6.  Nondimensionalized  mean  temperature  gradient  compared  with  LES  results. 
The  data  points  indicate  the  average  over  the  three  traverses;  brackets  indicate  the  stan¬ 
dard  deviation. 
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Figure  7.  Mean  dye  concentration  versus  height  at  five  times. 


Figure  8.  Time  rate  of  change  of  mean  dye  concentration  versus  height  at  five  times. 
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Figure  9.  Deduced  vertical  concentration  flux  profiles  at  five  times. 
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